Coulomb effects in semiconductor quantum dots 
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Coulomb correlations in the optical spectra of semiconductor quantum dots are investigated us- 
ing a fuU-diagonalization approach. The resulting multi-exciton spectra are discussed in terms of 
the symmetry of the involved states. Characteristic features of the spectra like the nearly equidis- 
tantly spaced s-shell emission lines and the approximately constant p-shell transition energies are 
explained using simplified Hamiltonians that are derived taking into account the relative importance 
of various interaction contributions. Comparisons with previous results in the literature and their 
interpretation are made. 
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I. INTRODUCTION 

Recent experimental investigations of optical spec- 
tra from individual self-assembled semiconductor quan- 
tum dots (QDs) [l| renewed the interest in funda- 
mental quantum-mechanical effects of confined interact- 
ingfew-carrier systems and their theoretical description 
Qi IE H> IE IS ' The few-carrier systems in individual QDs 
are interesting candidates for potential applications in 
quantum optics like sources of single photons or entan- 
gled photon pairs Q- 

Apart from the experimental relevance, the problem of 
few charge carriers in the discrete states of a given con- 
finement potential is also a paradigm for an interacting 
few-carrier system, that can be solved without further 
approximations. A discussion of the exact (but compli- 
cated) results in terms of simplified pictures is therefore 
desirable. 

In the past, electronic states and the resulting dipole 
transition have been calculated in box-like confinement 
potentials using exact diagonalization Q as well as in 
a two-dimensional system with harmonic confinement 
potential 3] with a limited number of configurations 
considered in the diagonalization. Furhermore, config- 
uration interaction calculations for a numerically deter- 
mined strain-induced confinement potential Q and den- 
sity functional calculations of the energy level structure 
and luminescence spectra [g] have been performed. 

In this paper we study the influence of Coulomb corre- 
lations on the optical spectra of QDs. Starting from the 
localized single-particle states in flat, cylindrically sym- 
metric QDs, described in the envelope function picture, 
the interacting states are obtained from the full diago- 
nalization of the Hamiltonian including Coulomb inter- 
action. From these states, optical emission spectra are 
computed for an increasing number of excitons assumed 
to be in their (interacting) ground states. Even though 
the Hamiltonian considered here is the same as in Ref . J3] , 
the full diagonalization results differ in several respects 
from the approximate picture obtained there. 

Despite the complexity of the spectra, characteristic 



features appear that can be explained in terms of simpli- 
fied Hamiltonians. These are obtained from the original 
Hamiltonian by retaining dominant terms and neglecting 
less important ones, such that, on the one hand, the es- 
sential spectral features are preserved and, on the other 
hand, analytic results can be deduced. The aim is to get a 
more intuitive picture of the influence of Coulomb effects 
and to provide an alternative to a numerically demand- 
ing full diagonalization approach. The trade-off between 
accuracy and simplicity can be reached in several ways 
and we give here two examples whose merits and short- 
comings are assessed against the full calculation result. 



A first case is represented by a diagonal Hamiltonian, 
which is obtained by retaining only the direct and ex- 
change terms, and which describes fairly well the main 
spectral features. The role of the Coulomb exchange in- 
teraction in the ladder-like structure of the lower part 
(s-lines) of the emission spectrum, recently observed in 
Refs. |j,l3 can be easily accounted for in this description. 



For explaining the absence of a similar ladder structure 
in the upper part (p-lines) of the spectrum a more careful 
handling of the Coulomb interaction is required. This 
is done in the second simplified approach described by 
an adiabatic Hamiltonian. One encounters here the so- 
called 'hidden symmetry'. This prop erty was discussed 
in a series of papers [IE [HI HE HEl- Nevertheless, the 
proof of the argument makes use of the assumption that 
one deals with a single degenerate shell, while in fact 
several, not weakly interacting shells are always present. 
Therefore it is by no means clear if the 'hidden symmetry' 
result holds, and if it does why and in what form. We 
show that in the adiabatic approach, and due to certain 
peculiarities of the problem, the 'hidden symmetry' can 
be recovered, albeit with the parameters renormalizcd by 
the interaction between the shells. The discussion of the 
conditions in which this happens sheds light on the limits 
of validity of the 'hidden symmetry' argument. 



II. HAMILTONIAN 

To describe the system of QD electrons and holes under 
the influence of Coulomb interaction we use the Hamil- 
tonian H = Hq + HqouI in the envelope-function approx- 
imation, 
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Here Cio- Wia) are annihilation (creation) operators of 
electrons with spin a in the one-particle orbital states |i) 
of energy e,^. The corresponding operators and single- 



and the superscript ee of the latter will be skipped in 
what follows. 

The Hamiltonian is fully determined by the effective 
masses me = 0.065 ttjq and rrih — O.ITtoq for the 
electrons and holes, respectively, the dielectric constant 
£,. = 13.69 as well as the oscillator length lose — 5.4 
nm (describing the in-plane confinement) and the QD 
height L = 4 nm. From these parameters we ob- 
tain the single-particle energies with constant spacing 
hwe = h^ /mellsc = 40.20 meV and hwh = 15.37 meV as 
well as the energy scale of the Coulomb matrix elements 
Ec — e^/ATreoCrlosc = 19.48 meV. In the following we will 
restrict ourselves to a QD containing only s- and p-shells 
for both electrons and holes. In this case the orbital of the 
single-particle state \i) can be uniquely identified by its 
angular-momentum, m = for the s-shell and m = ±1 
for the two p-states. 

The six-fold integral in Eq. |2Jl can be analytically 
reduced to a one-dimensional integral, which is deter- 
mined numerically. Due to the cylindrical symmetry of 
the problem we have angular momentum conservation: 
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particle energies for holes are h^„ (h]^) and £,'*, respec- elements are listed in Table U 
tively. The Coulomb interaction matrix elements are 
given by 



All non-zero Coulomb matrix 



V, 
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with the band index A — e,h, the single-particle wave 
function {r\i,X) — S,i^x{'r), and the Coulomb potential 
V{r) = e^ /AireQerr, where e^ is the background dielectric 
constant. 

It has been shown that in the case of flat, cylindri- 
cally symmetric QDs the single-particle bound-state wave 
functions in the plane of larger extension (perpendicular 
to the growth direction) are well approximated by those 
of a two-dimensional harmonic oscillator [lj| . Due to the 
rotational symmetry around the dot axis the correspond- 
ing angular momentum is a good quantum number. For 
the strong conflnement in the growth direction we use an 
infinite potential well to model the corresponding finite 
extension of the wave-function. Only the energetically 
lowest state due to confinement in the growth direction 
will be considered. 

A particular situation which is often employed in the 
literature and will be adopted in the following is the so- 
called symmetric case, in which one assumes identical 
envelopes for the valence- and conduction-band electrons, 
i-e. £,i,e = ^ih — ^i- This assumption holds exactly for 
QDs modelled by a box-like potential and is found to 
be a good approximation for oscillator potentials. Note 
that for the same state index i the electrons and holes 
have opposite angular momenta. With the assumption 
of identical envelopes all Coulomb matrix elements are 
related to the electron-electron ones 





{i,j,k,l) 




Vij^ki/Ec 




(0, 0, 0, 0) 




1.1197 




(0, 1, 1, 0), (1, 0, 0, 1) 




0.8690 




(0, -1, -1, 0), (-1, 0, 0, -1) 




0.8690 




(1, 1, 1, 1) 




0.7935 




(1,-1,-1, 1),(-1, 1, 1, -1) 




0.7935 




(-1, -1, -1, -1) 




0.7935 


(0, 0, 1, 


-1), (0, 0,-1, 1),(1,-1, 0,0), (-1, 


1, 0, 0) 


0.2507 




(0, 1, 0, 1), (1, 0, 1, 0) 




0.2507 




(0, -1, 0, -1),(-1, 0, -1, 0) 




0.2507 




(1, -1, 1, -1),(-1, 1,-1,1) 




0.1753 
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TABLE I: Non-zero Coulomb matrix elements in units of 
Ec- The indices refer only to the ^-component of the angular 
momentum m = for the s-shell and m = ±1 for the p-shell. 
The horizontal lines divide the matrix elements into three 
groups: direct- (top), exchange- (bottom) and other terms. 



III. EXACT DIAGONALIZATION 

In a semiconductor QD the finite height confinement 
potential leads to a finite number of localized states as 
well as to energetically higher delocalized states. When 
the influence of the delocalized states on the discrete QD 
spectrum is neglected, the eigenvalue problem has a fl- 
nitc (albeit large) dimension and can be solved without 
further approximations. In this way the Coulomb in- 
teraction between the different possible configurations of 
carriers in the available bound states is fully taken into 
account. 

As the Hamiltonian conserves the total number of elec- 
trons A^e and holes Nh, the Hamiltonian matrix falls into 



subblocks with basis states corresponding to uncorrelated 
many-particle states of the form 



n (4)"^ n 



{h]) 



|0) 
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Moreover, the Hamiltonian (^ commutes with the total 
(electron plus hole) angular momentum Iz =l% + l^ , with 
the total spin of electrons, S'^, S^, as well as with the total 
spin of holes, Sf^,S^. This rather rich symmetry can be 
used for separating even smaller Hamiltonian subblocks 
and for predicting degeneracies. A list of good quantum 
numbers includes therefore TVe, A^/i, Zz, S^, S^^, S'f , 5^ and 
the eigenstates are degenerate with respect to S^,S^. By 
numerical diagonalization of these subblocks one gets the 
Coulomb-correlated states and the corresponding eigen- 
values, classified according to the above-mentioned sym- 
metries. 

The coupling of the QD with the optical field is de- 
scribed by the interband dipole operator 



-P = ^ ^cr = rf ^ hi,aei-a 



(5) 



and its hermitian conjugate, where d is the interband 
dipole matrix element. The operator P describes the 
recombination of "mirror-symmetric" pairs, i.e., having 
opposite ^-components of the angular momentum and 
spin. The diagonality in the orbital index i is a conse- 
quence of the identical envelopes of the electron and hole 
states. The matrix elements of P define the QD emission 
as given by Fermi's golden rule. 






where the correlated initial (i) and final (/) states and 
their energies are calculated as described above. The final 
state has, of course, one electron-hole pair less than the 
initial one. The hermitian conjugate operator P^ creates 
mirror-symmetric pairs and appears in a similar formula 
for the light absorption. 

Next, one has to specify the initial states for which the 
optical spectra will be calculated. In the following we 
restrict ourselves to situations where optical excitation 
leads to the same number of electrons and holes within 
the QD, Ne = Nh = N^ where N^ stands for the num- 
ber of electron-hole pairs (in the following loosely called 
excitons). It is further assumed that energy relaxation 
of carriers within the QD is considerably faster than car- 
rier recombination such that the initial state for the re- 
combination process with given N-^ is the correspond- 
ing multi-exciton state with the lowest energy {N^ exci- 
ton ground state). Moreover, changes of the carrier spin 
during relaxation are neglected [lj|. Since optical exci- 
tation only creates multi-exciton states where electrons 
and holes have opposite S^ values, only states with van- 
ishing total z-component of the spin are considered as 
initial states in Eq. ®. 



The analysis of Coulomb-correlated multi-exciton 
states shows that for even N^ , the symmetry of the 
ground state is singlet-singlet (ss), i.e., Se = Sh = 0. 
These states are nondegenerate. On the other hand, for 
odd N-^ one has doublet-doublet (dd) ground states with 
Se = Sh — 1/2 which are four times degenerate. For the 
above discussed choice of the initial states for the recom- 
bination process only two states contribute to the emis- 
sion formula with a weighting factor of 1/2 each. The 
other two states are also dipole-allowed, but are elimi- 
nated from the emission formula because their total Sz 
is nonzero. (By the same arguments one may be con- 
cerned about eliminating the states not having total an- 
gular momentum zero but, as expected, no such ground 
states occur.) The final state can be any (ground or ex- 
cited) state of the system with one exciton removed. 
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FIG. 1: Ground state emission spectra for a quantum dot 
with different number of excitons. The labels indicate the to- 
tal spin of the final state for electrons and holes, respectively: 
s-singlet, d-doublet, t-triplet. The total angular momentum 
z-component is always zero. All energies are measured rela- 
tive to the band-gap. 



The results of an evaluation of Eq. © based on Cou- 
lomb-correlated multi-exciton states are shown in Fig. ^ 
for an initial filling of one up to six electron-hole pairs. 
For better visibility the J-functions are represented by 



narrow Lorentzians centered at the energy hiu — Ei — Ef 
(measured relative to the band gap) and having the area 
equal to the oscillator strength |((;z!i/|P|(^i)p . 

Several features are immediately obvious. First, there 
is a clear spectral separation between the higher energy 
p-lines for fku > 80 meV (obtained by removing one ex- 
citon from the p-shell) and the lower energy s-lines (for 
fko < 35 meV). Second, the position of the s-lines are ar- 
ranged approximately in a descending ladder as the num- 
ber of excitons increases, while the energy of the p-lines 
show a remarkable stability. The last fact was attributed 
to a 'hidden symmetry' property, which will be discussed 
below. The dipole operator P has no simple commuta- 
tion relation with the spin operators Sg, S^^ and therefore 
the spin symmetry of the final state is not determined by 
that of the initial state [16] . This is why one encounters 
as final states all possible spin symmetries (ss, dd, tt, st, 
ts), as indicated in the figure. 

By restricting the present discussion to s- and p-shells 
only, the Hamiltonian of Eq. ^ is identical (up to 
nonessential differences in the parameters) to the case 
analyzed in Q. Nevertheless, the full diagonalization 
procedure used here leads to different relative line inten- 
sities and, more importantly, to the appearance of new 
emission lines (the st- and ts-lines are missing in [3). 



IV. APPROXIMATE HAMILTONIANS 

The interaction between different state configurations, 
as given by the exact diagonalization procedure, shows a 
quite high degree of complexity and therefore the results 
are not immediately intuitive. Even though we have con- 
sidered electron-hole pairs that are optically created in 
mirror-symmetric states, the Coulomb interaction mixes 
them strongly with configurations in which the electrons 
and holes are not arranged symmetrically. For instance, 
promoting two holes from the s-shell to the p-shell is en- 
ergetically less costly than promoting one electron and 
one hole. The second case is more symmetric, but the 
first produces a state which is energetically closer and 
therefore participating stronger in the exact interacting 
eigenstate. This may explain the disagreement with the 
line intesities found in 'j||, where only symmetric states 
are considered. 

On the other hand, the relatively regular structure of 
the emission lines seems to indicate that an intuitive pic- 
ture should be possible. This is achieved by turning to 
approximate, simpler Hamiltonians which allow analytic 
solutions and at the same time retain the essential fea- 
tures of the full problem. 

One such simplified Hamiltonian can be obtained as 
follows. By examining Tabled one sees that the largest 
Coulomb matrix elements are the direct ones, Vij,ji = 
Dij = Dji. Their contribution to Hcoui can be ex- 
pressed solely in terms of the occupation number opera- 
tors, nf^ = ej^ejo- and fi'^^ = hj^hi^, and therefore is di- 
agonal in the noninteracting basis, Eq. Q. The same is 



true for the exchange matrix elements VijAj = Xij = Xji 
with i ^ j provided one includes their contribution only 
in the e-e and h-h interaction terms involving the same 
spin {a — a'). In this way one obtains the diagonal 
Hamiltonian 



H, = 



1 






-Du I n. 
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l,3,cr 



The prime in the last summation indicates that i = j 
terms have to be omitted and n^' = n|j -(- n^\ . Of 
course, for this Hamiltonian there is no configuration in- 
teraction. The non-correlated states are eigenstates and 
the eigenvalues are derived by the above formula by in- 
serting the corresponding occupation numbers. 

If the states are symmetrically populated, n^^ = 



.X 



their energy is further simplified: 
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D,,. 
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(8) 



In this model, the exciton energy is Ef = ef + e^ — Da 
where the binding energy results from the direct electron- 
hole Coulomb attraction. When the electron and hole en- 
velope functions are identical, the excitons are not only 
globally but also locally neutral and the direct electro- 
static interaction between different excitons vanishes. In 
the approximate Hamiltonian ^ the only interaction be- 
tween excitons comes from the exchange terms and takes 
place between excitons with the same spin structure 17]. 
For instance, the biexcitonic binding energy is zero in 
this approximation (the full model gives a small binding 
energy of about 2 meV). A comparison of the diagonal 
model with the full result is given in Fig. |21 The ap- 
proximate spectrum indeed shows a ladder-like structure 
of the s-lines, with the corresponding spacing in good 
agreement with the full calculation and given by the ex- 
change interaction between carriers in the s- and p-shell, 
2Xsp- According to Eq. ((SJ the energy change for the re- 
moval of one s-exciton with a given spin structure is pro- 
portional to the number of p-excitons having the same 
spin structure. Adding equal contributions from elec- 
trons and holes the coefficient is 2Xsp. Also, in the case 
of an odd number of excitons one has two s-lines, depend- 
ing on whether the spins in the removed s-exciton agree 
or not with the spin orientation in the majority of the 
p-excitons. For an even number of excitons this splitting 
does not occur. These main spectral features have been 
described previously |j, |2| , here we show which terms of 
the full Hamiltonian are responsible for them and that 
these terms can be included in an exactly solvable ap- 
proximate Hamiltonian H^. 

Even though the diagonal Hamiltonian gives an intu- 
itive picture of the main features in terms of uncorrelated 
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FIG. 2; Comparison of results from the diagonal Hamilto- 
nian, Eq. Q, (solid line) and the full diagonalization (dashed 
line) for quantum dot emission spectra with increasing num- 
ber of excitons. 



states, the following shortcomings of this model need to 
be discussed. The ground state of the 4-exciton prob- 
lem is inaccurately given by the diagonal model as a 
triplet-triplet state, as predicted by Hund's rule. It is 
known that Hund's rule does not always apply in QDs 
and this is such a case. The true ground state with 
an energy slightly below this triplet-triplet state has a 
singlet-singlet symmetry, as mentioned in the previous 
section. It was this singlet-singlet state which was used 
as initial state in the emission spectrum obtained from 
the diagonal Hamiltonian, solid line in Fig. [3 in order 
to make the comparison meaningful. Also the diagonal 
Hamiltonian entails a ladder-like structure for the p-lines 
too, with a spacing of 2Xpp, which is not confirmed by 
the full calculation. As discussed below, this stems from 
neglecting important pieces from the interaction. 

Regarding the p-shell emission, one can refine the ap- 
proximate Hamiltonian in the following way. An exam- 
ination of the ground states for a situation with three 
or more excitons, calculated with the full diagonaliza- 
tion procedure of Section IhH shows that it is safe to as- 
sume that in such states the s-shell is completely filled. 



The states with full s-shell configurations appear with 
a weighting factor of at least 0.95. Therefore, as far as 
the ground states and the lower excited states are con- 
cerned, it is possible to construct an approximate Hamil- 
tonian describing the fully interacting p-states following 
'adiabatically' an external field provided by the 'frozen' 
s-carriers. Practically this is obtained from Eq. Q along 
the same lines as before, but this time one enforces di- 
agonality only with respect to the s-occupation numbers. 
In other words, one discards only those terms which con- 
tain s-state creation or annihilation operators and cannot 
be expressed in terms of s-state occupation numbers. In 
this way, one obtains fully correlated p-states at given 
s-orbital fillings. 

In the resulting 'adiabatic' Hamiltonian Had = Had + 

■^ad + ^ad ^^ havc Separated the terms describing the 
s- and p-shell as well as the s-p-interaction. The p-shell 
part retains the form of Eq. ^ with the summation re- 
stricted to the p-orbitals. Therefore, in the following 
equations, the indices i,j,k,l label only p-states while 
for the s-states the explicit subscript s is taken. Using 
the symmetry relations of the Coulomb matrix elements, 
Eq. Q, and bringing close the operators with the same 

(v) 
spin, H^J can be rewritten as 

h':J = E^?"^-^E^.-^.4^ 

i(T ijka 
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X (e]_^,efe,^'-4_^,%„.) (9) 

where the new one-particle Coulomb terms result from 
the reordering of operators. The s-shell part and the 
s-p-interaction are similar to the diagonal Hamiltonian, 

Eq. 0, 



^efco 



^'y.fcj ^lo-Zifccr 



a 



ad 



-D, 



2 



2 " 



n^^ 



H. 



(sp) 



D, 



(-; 



X, 



■E( 



^s,CT^p,cr 



^s,(j^^p,a) 



(10) 



(11) 



where n^'^ stands for the total population of electrons or 
holes on the s- or p-shell, respectively. It is important 
to note that, with the Coulomb matrix elements listed in 
TableU in the one-particle Coulomb terms of Eq. (|5J only 
the contributions with i = k appear. They are responsi- 
ble for a renormalization of the one-particle energies e^' . 



Moreover, these renormalized energies do not depend on 

(v) 
the index i so that H^J describes a 'single degenerate 

shell' l^ilijjll^, and can be rewritten as 
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With the help of the adiabatic Hamiltonian Had the 
full emission spectrum can be explained as follows. For 
the s-lines the arguments showing the formation of a lad- 
der with the spacing of 2Xsp are as in the case of the 
diagonal Hamiltonian Hd- The energetic position of the 
p-shell emission can be deduced from the commutation 
relation of Had with the p-shell dipole-transition opera- 
tor Pp, defined as in Eq. JSJ but with the sum restricted 
only to p-states. It is readily verified, that 



A 



h^ h P 







(13) 



and therefore Pp commutes with the last two lines of 



Eq. p2ll . This is the core of the 'hidden symmetry' prop- 
erty 10,, 13j showing that the interaction part plays no 
role in this argument. The commutation of Pp with the 
occupation number operators of the shells is rather ob- 
vious, because the application of Pp corresponds to a 
population reduction by one electron-hole pair in the p- 
states and no change in the s-states. Correspondingly, 
one can verify [hp'^, Pp] = ~Pp and [hf'^, Pp] = 0. Such 
simple relations arise only when commuting Pp with the 
total number operator of the p-shell, not with individual 
number operators. This is why it is important to have 
degenerate shells. 

Assuming, that the s-states are fully occupied 
{nf'^ = 1) one obtains from these results 



[Had, Pp] 



\^p I ^p ^pp ^pp 
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(14) 

This shows that the removal of one p-exciton is accom- 
panied by an energy decrease which does not depend on 
the number of excitons. The value of this energy, 
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is in excellent agreement with the exact diagonalization 
result. Indeed, using the values in Tabled one obtains 
£p = 82.5 meV. The emission spectrum associated with 
the adiabatic Hamiltonian consists therefore of p-lines 
having all this common value, while the s-line ladder re- 
mains the same as given by the diagonal Hamiltonian. 



The 'hidden symmetry' argument, as discussed by 
Wojs and Hawrylak, 0,113, lldl is proven on the assump- 
tions that (i) one has identical envelopes in the two bands 
and (ii) the one-particle levels form 'a single degenerate 
shell'. In these hypotheses one gets a constant energy 
value at the removal of each exciton, and this value de- 
pends only on the parameters of this shell. 

Nevertheless, the Hamiltonian under discussion here, 
as generally used in the literature, is rather describing 
several interacting shells, so that is not obvious why (if 
at all) the argument holds. In the present case the answer 
is contained in the adiabatic Hamiltonian. In it the p- 
shell is the 'single degenerate shell' because (i) the s-shell 
is 'frozen' and higher shell are absent in the considered 
situation and (ii) the field created by the s-shell carriers 
does not remove the degeneracy of the p-shell. In such 
conditions the adiabatic Hamiltonian obeys exactly the 
'hidden symmetry' commutation relations. In this pic- 
ture the energy for the removal of a p-exciton, Eq. H15|l . 
contains also terms coming from the s-p-interaction, and 
it is this value that is in agreement with the full diago- 
nalization. 

This discussion also shows the validity range of the 
'hidden symmetry'. For instance, the presence of higher 
shells (but also depending on the actual model parame- 
ters) may spoil the argument. This seems to be the case 
described in |j|, where the p-lines are not independent 
on the exciton number, but are arranged in a descending 
ladder too. Also, the field of the 'frozen' states may re- 
move some of the degeneracy of the shell in question. For 
example if the outer shell is a d-shell, the states with zero 
angular momentum and those with angular momentum 
±2 will experience the field created by the s-shell carriers 
differently . 

One may argue that by approximating the full Hamil- 
tonian, Eq. (Q, with the adiabatic one, Eqs. l|^- Hlll) . the 
interaction inside the p-shell is still too complicated to al- 
low analytic diagonalization, i.e. the 'hidden symmetry' 
property is a simple relation between otherwise complex, 
strongly correlated states. Nevertheless, the Fock sub- 
space generated by the p-orbitals is significantly smaller 
and this is in itself a numerical simplification. Symmetry 
arguments also can be used to reduce the blocks to be 
diagonalized and at least the ground states can be ob- 
tained analytically. A procedure for obtaining analytical 
eigenstates is the repeated applications of the raising op- 
erator Pt o n the 'vacuum' (full s-shell, empty p-shell) 
state HESIll- 

In summary, for semiconductor QDs with finite height 
two-dimensional harmonic confinement potential the 
multi-exciton emission spectra are discussed on the ba- 
sis of a full diagonalization of the Hamiltonian including 
Coulomb interaction for the localized states. The charac- 
teristic features are s-shell lines arranged approximately 
in a descending ladder with increasing exciton number as 
well as nearly constant energetic position of the p-shell 
lines, provided that these are the only confined shells. 
Based on the relative importance of the various Coulomb 



matrix elements, a simplified Hamiltonian has been con- 
structed which is diagonal in the single-particle states. 
It explains the s-shell emission properties, that appear 
as soon as the p-shell population starts to contribute, in 
terms of the s-p-exchange interaction Xgp. The results 
are two main s-shell lines for odd number of excitons 
separated by 2Xsp, while the spectrum is dominated by 
a single line for an even exciton number. The almost con- 
stant energetic position of the p-shell emission ('hidden 
symmetry') is discussed in terms of Coulomb correlated 
p-shell carriers in the presence of a completely filled s- 
shell. It is also shown that the arguments for the 'hidden 
symmetry' break down as soon as higher confined shells 
contribute. Depending on the coupling strength of higher 
shells to the p-shell the energetic stability of the p-lines 
is expected to disappear gradually. 

The proposed simplified Hamiltonians give a more in- 
tuitive picture of the rich properties of the emission spec- 



tra from the Coulomb-correlated QD carriers. They also 
might be an alternative to a demanding full diagonaliza- 
tion scheme when Coulomb interaction in the presence of 
three and more electron-hole pairs is important. Finally, 
only a diagonal Hamiltonian allows to discuss physical 
processes in terms of single-particle states. Our compar- 
ison of results from full and diagonal Hamiltonian shows 
to what extent this is justified. 
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